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Abstract 

We consider a homogenization problem for magnetic Ginzburg- 
Landau functional in domains with large number of small holes. For 
sufficiently strong magnetic field, a large number of vortices is formed 
and they are pinned by the holes. We establish a scaling relation 
between sizes of holes and the magnitude of the external magnetic 
field when pinned vortices are multiple and their homogenized density 
is described by a hierarchy of variational problems. This stands in 
sharp contrast with homogeneous superconductors, where all vortices 
are known to be simple. The proof is based on F-convergence approach 
which is applied to a coupled continuum/discrete variational problem: 
continuum in the induced magnetic field and discrete in the unknown 
finite (quantized) values of multiplicity of vortices pinned by holes. 

1 Introduction 

Vortices determine electromagnetic properties of superconductors that are 
important for practical applications (e.g., resistance). A key practical issue 
is to decrease the energy dissipation in superconductors, which occurs due 
to the motion of vortices. This dissipation can be suppressed by pinning 
of vortices. In particular, a physical problem of pinning in superconducting 
thin films with a periodic array of antidots (holes) was considered in [?] (see 
also references therein) . This problem leads to analysis of a two-dimensional 
Ginzburg-Landau (GL) energy functional in a domain with periodic holes. 

In this work we consider a mathematical model of pinning of vortices by 
many holes in relatively small superconducting samples (comparable to the 
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London depth) submitted to a uniform magnetic field, which is weak so that 
the vortices do not appear in the bulk of the superconducting sample and 
they may exist only in the holes. 

Since modern experimental techniques allow to create very small holes, 
the question arises what kind of physical effects can be expected for such 
sample. In particular, typical experimental results lead to uniform arrays of 
vortices in the entire domain. In this work we present a mathematical model 
that leads to special "critical" scaling when a nested sequence of sub domains 
with vortices of increasing multiplicity appears. 

We next present a brief review of relevant mathematical work. The study 
of pinning by a finite number of pinning sites was pioneered [10] where a 
simplified GL model (SGL) with no magnetic field was and discontinuous 
pinning term for a single inclusion was considered. The existence of d vorticies 
of degree 1 inside the inclusion was established when Dirichlet boundary 
data with degree d. The results of [10] were subsequently generalized for 
the magnetic GL functional [8], [9] and pinning by a single inclusion. A 
comprehensive study of pinning by finitely many normal inclusions and holes 
for the the magnetic GL functional was performed in [H [2] . More recently 
pinning by finitely many holes whose sizes goes to zero as the Ginzburg- 
Landau parameter goes to infinity was established in [5] for the SGL model. 

Homogenization in the framework of magnetic GL model with continuous 
oscillating pinning term was considered in the pioneering work [3] , where large 
number of vortices is described by the homogenized vorticity density. Since 
some composite superconductors are described by a discontinuous pinning 
term, in subsequent works [6l [7] homogenization problems for such term were 
address in the context of simplified GL model and special Dirichlet boundary 
conditions, which result in either no vortices [B] or c/ vortices [7] . 

In this work we study a homogenization problem for a large number of 
vortices and large number of pinning holes that are described by a perfo- 
rated domain (which corresponds to a discontinuous pinning term). This 
problem is described by the minimizers of the GL functional 

GL[u,A] = l [ {\\/u-iAu\^ + ^{l-\u\'^f)dx+l- f {cm\A-hUdx. (1) 

The unknowns here are the complex order parameter u and the vector po- 
tential of the magnetic field A, while h^^^ is given external magnetic field 
(positive scalar number). The domain fi^ in ([1]) is obtained by perforating 
a given simply connected bounded domain f2 C by a large number 
of small holes. Holes are all identical disks with radius and periodically 
distributed centers a^, with e > being a small spatial period. We assume 
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that 

I logft;| > I log Pel and Pe <ti e (2) 

(holes radius much greater than vortex core and much less than the spatial 
period e). Moreover, we consider the following scales of magnetic field and 
diameters of holes, 

he.t = \/e\ diam(a;j) = 2e~'</'\ (3) 

where A and 7 are fixed positive numbers. 

Formal calculations show that under the scale separation condition ([2]) 
there are no vortices in the bulk of the domain Vt^ and the term with the 
integrand ^(1 — \u\^Y can be effectively replaced by the constraint \u\ = 1. 
Thus we are led to the minimization problem 

M, = M{Fe[u, A]; ue H^Vle; S^), A G H\n; R^)} (4) 
for the functional 

F^[u, A] = - I I Vu - iAu\^dx + - / (curM - K^^fdx, e > 0. (5) 

For every fixed e > minimizers of (m^,/!^) of problem (ll])-© can be ex- 
pressed in terms of degrees of the order parameter on the boundaries of holes 
as follows, see Let c?^ be (integer) degrees of on (9a;j, then the induced 
magnetic field = curM^ satisfies 

' -Ah' + = 0m n„ 
h'{x) = /ig^j on d^, 
}f{x) = H^ on cui, j = l,2,...,iV„ 



-L^/£d^ = ^^d^-L^Hx)dx, 

where Hj are unknown constants that are part of the problem. Note that if 
we know the degrees dp then the minimum (jlj) is given by 

Esihe) = \l \Wh'\'dx + \l {h' - hWdx, (7) 

and conversely, the minimum of (E]) is obtained by minimizing ([7]) in integer 
dp where /i^(x) defined via {d'j} as the unique solution of (jH]). 

Our goal is to describe the asymptotic behavior of the degrees as e — )■ 0. 
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2 Homogenization (corrector) and compact- 
ness results 



We introduce rescaled quantities 

= e^h\ = e^hl^t, EeChe) = e^E,(he), M, = e^M,. (8) 

Note that (and therefore /i^) is determined uniquely by the tuple of inte- 
gers (ij. Thus, abusing a little notations, we may write Eeijf) = E^dd'^j}). 
Consider a minimizing tuple of degrees {(i^}, so that the corresponding solu- 
tion of Qj, rescaled according to ([8]), satisfies Ee{hs) = M^. 

First we obtain a priori bounds for these degrees dj in the following 

Lemma 1. Let dj be degrees of the minimizer of then 

J2id^r < Cle\ (9) 

where C is independent of e. 

Proof. The weak formulation of the problem for reads, find G H^{VL) 
such that V/i^ = in all cjj and /i^ = A on and 



I {Vh' ■ Vv + h'v)dx -e^Y^ d)v\^e = 



(10) 



holds for every test function v G Hq{Q) such that Vf = in all Uj. In 
particular, if we choose all dj = 0, set v = }f — \, we get the a priory bound 
(for h corresponding to this choice of degrees dp 

:= [ {\Vh'\^ + Ch'y)dx < c 

therefore < C, where C is independent of e. Hence for the minimiz- 
ing tuple {c?^} we have -Ee[/i^] < C and thus ||/i^||_f/i(f7) < C, with an- 
other independent of e constant C . Now choose the test function v = 
Y,d)L){x)/\og{2pe/e) in where 



^ ^ \og{2pJe) if X G 5p^(ap, and if x ^ 5,/2(ap. 
Then simple computations lead to the required bound, 

2 ^(dp2 = I (y^e . y^, ^ j^e^,^^^ < ^2716^1 + o{e)) Y.id'^VlWWm 

i.e. £2E(rfp'<27r(l + o(5))||/.^||^,(^)/7<a □ 
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It follows from Lemma [T] that, up to extracting a subsequence, 

C = e^^ d]5a-{x) D{x) as distributions, and D{x) E L'^iVt). (12) 

From now on {d'f\ denotes an arbitrary sequence of tuples of integers such 
that ([9]) and f lT2|) hold, and h is the function associated to the tuple {d!j}, i. 
e. /i^ = e^/i^, where /i^ is the solution of 

It is rather easy to see that under the above mentioned conditions we can 
pass to the limit in f IlOp to get that converges if^-weakly to the solution 
h of the homogenized problem 

-Mi + h = 2TiDix) mVt , , 

(13) 

/i = A on d^. 

(for details see Lemma [3] below). However this result is too week for our 
principal goal of describing the limiting vorticity D{x), which will be (natu- 
rally) done by calculating the F-limit of the functionals -E^. In the next step 
we introduce a corrector and get the strong if^-convergence. 
We consider the ansatz, 

h%x) = h%x) -e^Yl d^j^ix) = h%x) + i^^ (14) 

where functions L'^j{x) are given by flllj) . The problem for /i^ (in its weak 
form) is, find Jf G H^{VL) such that Vh^ = in all and Jf = \ on dQ, and 



[ {Vh' ■ Vv + hev)dx + V / vR 



[x)dx 



Y,d] [ vds (15) 



holds for every test function v G Hq{Q) such that Vv = in all Uj. 

Lemma 2. Under condition ( f7^) functions converge H^-strongly to h, the 
unique solution of [Tl 



Remark 1. Lemma\^ shows that the function R^{x) = —e'^'^djLj{x) is a 
corrector, so that h'^ = — converges strongly in H^{Q). 

Proof. Using Lemma [1] one shows that i?^ converges if ^-weakly to zero, and 
therefore h'^ h, up to extracting a subsequence . To prove that h solves 
(fT3|l we consider test functions v'^ = v{x) + ^ 0(x/pe)(u(ap — v{x)), where 
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V G C^{VL) is an arbitrary function, and </'(x) is smooth cut-off function such 
that = 1 if |x| < 1 and = if |x| > 2. Set v = in f[T^ and pass to the 
hmit as e —7- to get 

{Vh-Vv + hv)dx = 2Ti / D{x)vdx. 



Thus h solves f[T3|) . 

Next we show that /i^ converges iJ^-strongly to h. We set w = /i^ — A to 
obtain in the hmit e — )■ 0, 



hmsup / V/i^ ■ V/i^(ix = hmsup5 yjc/^ / {h^ — X)ds+ I {X — h)hdx. 

By the Pioncare inequahty, we have 

/ h'ds = \ [ h'dx + 0{l){[ \Vh'\^dxY^\ 

where 11^ is the cell with the center at and the side length e. Therefore 
limeVjc?^ / If ds = 2ti I D{x)hdx, 

where we have used Lemma [H (IT^ and the fact that h'^ ^ h strongly in 
L'^{Q). Thus, taking into account fll3p . we finally get 

limsup / W h"" -y If dx = 2tx i D{x){h-\) dx+ I {X-h)hdx = / Wh-Whdx. 
Jn Jn Jn Jn 

This implies that h'^ ^ h strongly in H^{Q). □ 

As a corollary of Lemma |2] we obtain 
Lemma 3. Under condition [W^) the following energy expansion holds, 

E,{h') = E,{h) + T^ie' + o(l), (16) 

where 

Ei{h) = l [ \Vh\^dx + l [ {h-\fdx (17) 



Proof. Since h'^ ^ h strongly in H^{VL) while i?"^ — ?■ weakly in if^(f2), we 
have 

Ee{h') = EiQi') + ]- I \VR'\'^dx + o{l) = Ei{h) + \ I \VR'\^ dx + o{l) 
^ Jn ^ Jn 

The straightforward calculation of the second term in this expansion yields 
(US]). □ 
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3 Limiting vorticity via F-convergence 

The main result of this work describing the hmiting vorticity is obtained by 
proving the F-convergence of functionals with respect to weak convergence 
f[T2|) of vorticity measures, 

E, r-converge to Eo{D) = Ei{h) + n-f (S?{D{x))dx as e ^ 0, (18) 

Jn 

where h is the unique solution of (fT3|) . More precisely we demonstrate that 

(i) (r — liminf inequlity) if conditions and f[T^ are satisfied then 

InnM E,{{d'^}) > Eo{D); (19) 

(ii) (r — limsup inequlity) \fD{x) G L^{Q) there is a (recovery) sequence 
of tuples {dj} satisfying conditions and (IT^ and such that 

limsup ^,(^^}) <Eo(D). (20) 

The function $ in the limit functional f[T5]) is a continuous piecewise linear 
function such that $(o?) = d"^ at integer points d. It describes the homoge- 
nized density of energies of individual vortices, whereas Ei corresponds to the 
interaction a vortex with magnetic field due to other vortices and external 
field. 

Thanks to the energy expansion 0161) . for "liminf" inequality we need 
only to prove the lower bound 

liminf£2^(d^)2 > 7r7 j $(D(x)) rfx. (21) 

Since the left hand side of (l^Tj) is nonlinear (quadratic) function of dp we 
use an analog of Young measures (see Appendix). 

3.1 Lower bound 

Spread the measure C^, which is the sum of point masses, over periodicity 
cells by setting = in 11^ (11^ is the cell centered at ap. Then 
represent as 



D%x) = J2 kf^li^)^ where /i|(x) = \ ^ ";7' (22) 
~Z 0, otherwise. 



1 in m, if A; = rf! 
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We extend and /c 7^ 0, on by setting = = in \ UII^ 
and also set fiQ = 1 in Q \ Ull^. The functions satisfy > and 

^ /i| = 1 and therefore form a partition of unity. Clearly, we can extract a 
subsequence such that 

^ /ifc weakly in L^(fi) VA; G Z. (23) 
Thanks to the bound (19]) we have 



Hence the limit functions Hk also form a partition of unity, 

/ifc > and /^fc = 1- (24) 
Moreover, the function D[x) defined in f|T2|) admits the representation 

D{x) = Y,kf^k{x) (25) 
and it is easy to see that the following inequality holds, 

liminfe^ J](ci^=)2 j /ifc(x). (26) 

In order to obtain a lower bound in terms of D{x) we prove the following 
simple 

Lemma 4. Given D (real number), then 

$(D) = min J ^ /cVfc; /ifc > 0, ^ /ifc = 1, ^ /c/ifc = D i (27) 
Ifcez J 

$(D) = (2A; + 1)|D| -A;-P ^//t < |D| < + A; = 0,1,2,. . . . Moreover, 
if D = d is an integer then then the unique minimizing tuple is = 1, 
rrik = ^ d. In the case D is non integer, represent D as the convex 
hull of two nearest integers d and d + 1, D = ad + {1 — a){d + 1), then the 
unique minimizing tuple is fid = a, fid+i = 1 — a, /^fc = Vfc ^ {d, d + 1}. 

Proof, li D = d is an integer then, clearly, $(-D) < d"^, and by Cauchy- 
Schwarz inequality we have 

d' = (J2kfi,Y <J2k'fik = HD). 



Thus $(-D) = (P and the minimizing tuple is = 1, = "^k ^ d. 

Consider now the case when D is non integer and D = ad+{l — a){d + l). 
Let {/ifc} be a minimizing tuple. If /ifc > for some k < d then there 
is jj^i > for some / > d + 1. Decrease fi^ and fii by a sufficiently small 
5 > and increase Hk+i and /i^-i by 6. This modification change neither 
nor (12^ but decrease the value of functional ^ k'^^k- Therefore /ifc = 
VA; ^ {d, d + 1}. The case when fi^ > for some k > d + 1 is similar. It 
follows from fl2^ and fl25l) that fid = C(, fJ'd+i = 1 — « and straightforward 
calculations yield the result. □ 

3.2 Upper bound 

In order to complete the proof of F-convergence (|T8l) we have to show the 
lim sup- inequality, i.e., given D G L'^{fi), we need to construct a sequence of 
tuples {d'^j} satisfying the boundedness condition IQ, that converge to D{x) 
in the sense of (fT2l) and satisfy inequality (!20|) . 

The limiting functional Eq{D) is continuous with respect to the strong 
convergence in therefore it is sufficient to establish the lim sup-inequality 

for D G C^{n) and then use density of C^i^) in L'^{n). 

Note that we not only need the converge of tuples in the sense of f|T2|) but 
more importantly we need the convergence of energies which does not follow 
from f[T^ . The key issue in the construction of the upper bound is that 
different configurations of vortices may lead tho the same vorticity D{x) = 
^kfik{x), however these configurations can be distinguished by ^k'^Hk{x) 
(which is equal to $(-D) for optimal Hk given by Lemma Hj). Thus we need 
to chose dj that define /i| via (122|) so that the limiting values fik are optimal 
in the sense of (12 7p . More precisely, if we represent D{x) for fixed x G as a 
convex hull of the nearest integers d and d + 1, D{x) = ad + {1 — a){d + 1) , 
we must have only holes with degree d and d+1 in a small neighborhood of x 
and in this neighborhood #{holes with degree d} / 7^{holes with degree 
should be approximately equal to a/ (a + 1). 

Recall that the centers of holes form e-lattice and therefore partition 
the domain Q into fine squares 11^ of side length e. We choose fixed sufficiently 
large integer M and consider squares Kk with side length (2M + l)e and 
centers at points a| which form an (2M + l)e-periodic lattice. Consider a 
square Kk strictly included in Q (it contains exactly (2A^ + 1)^ holes). Let 
Dk denote the mean value of D{x) on Kk, 



and let dk be the integer such that dk < Dk < dk + 1- Represent Dk as the 




9 



convex combination of (ifc and dk+1, = akdk+{l—<yk){dk+^) (0 < Ok < !)• 
Next we chose the largest integer R > such that R/{2M + 1)^ < Ok and 
set d'j := dk for R holes Uj in Kk and d'j := rf^ + 1 for other holes in Kk- 
We repeat this procedure for all squares Kk lying strictly inside fl and set 
degrees of remaining holes to be zero. Let functions be defined by (l22l) . 
then we have 

We claim that for sufficiently small e > 

E,Ch'))<Eo{D) + 6M, (28) 

where — >■ as M — )■ oo. Indeed, up to extracting a subsequence, fif — )■ 
///(x) weakly in L^(f2) for all / G Z. The limiting vorticity D(x) (which 
depends on M) is given by D{x) = '^lfii{x). Due to the construction of 
tuples {d'^j} we have, |/Xd(x) — a\ < 1/(2M + 1)^ and |/id+i(x) — (1 — a)| < 
1/(2M + 1)2 for every point x E n such that dist(L'(x), Z) > 1/(2M)2, 
where ad+{l — a){d+ 1) = D{x) is the representation of D{x) as the convex 
combination of nearest integers. On the other hand |/id(a;) — 1| < 2/{2M)'^ 
ii X E Q and \D{x) — d\ < 1/{2M)'^ for some integer d. It follows that 

lim£2^(dp2 = /" ^^(2;) < f ^D{x))dx + C/M\ 

We also have \D{x) — D{x)\ < C/M'^ in Vt. Thus, by virtue of Lemma 
[3] we obtain (128|) . Next choosing a suitable sequence of increasing integers 
M = Mir we get the required upper bound. 

3.3 F-convergence theorem 

We summarize results of this Section in the following 

Theorem 1. The junctionals Ef^ T-converge to Eq{D) as £ — )■ 0. The limit 
functional Eq{D) is given by 

Eo{D) = l [ {\Vh\^ + {h-Xy)dx + TT^ [ ^D{x))dx, (29) 
2 Jn Jn 

where h = h{D) is the unique solution of l[T3\) and $(-D) = {2k+l)\D\—k — k'^ 
if k < \D\ < k + 1, k = 0,1,2, ... , see FzgUi 

This yields the main homogenization result of this work 
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Figure 1: Function $(-D) 



Theorem 2 (homogenizad vorticity). Let {dj} be a tuple of integer degrees 
solving the minimization problem ^ (i.e. the solution of problem ^ 
minimizes the functional ^ ), then 

(x) ^ -D(x), in the sense of distributions, (30) 

where D is the unique minimizer of the functional Eq{D) in Eq{D) 
being given by (E^j- 

Proof. Note that the F-hmit functional Eq{D) is strictly convex continuous 
and coercive, therefore it has the unique minimizer. On the other hand, by 
Lemma [1] the weak limit e'^Y^^j^aji^) ~^ D{x) exists (up to extracting a 
subsequence) and D G L^(f2). Therefore, due to the classical properties of 
F-convergence, D is the unique minimizer of Eq{D). □ 
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4 Analysis of the limit problem via convex 
duality. Hierarchy of multiplicities 

We use convex duality (see, e.g., [12]) to pass from problem 

MA=min|^ / {\V{h-X)\^ + {h-Xf + 2TT^m-Ah + h)/{27i)))dx; 

{h- \) e H^iQ), -Ah + he L\Q)^ (31) 

to the dual one 

-M, = min|i^(|V/p + /2)dx + ^*(-/); /Gi/o'(^)}, (32) 
where J^*{f) is the Legendre transform of the functional 
J^(ft) = 7r7 / $((-Aft + fi: + A)/(27r)), 

i.e. 

J^*{f) = sup jy (V/ ■ Vk + fn) dx + J-(k); k E H^{n)^ . 

Due to the fact that J^{k) is lower semicontinuous, the minimizer (h — A) of 
(13T|) and minimizer / of (132|) coincide (moreover Mx in (l3Ti) and (l32l) is the 
same). Thus, we have 

27rD{x) = -Ah + h = -A/ + / + A (33) 

The calculation of the Legendre transform J^*{f) is reduced to the calcu- 
lation of the Legendre transform $* of the function 7r7$(K/(27r)). Indeed, if 
we use integration by parts we derive 

/ (V/ ■ Vk + fn) dx + T{k) = 
Jn 

[ {-Ak + K + X)fdx+ I {Tl-i<^{{-AK + K + \)/{2Tl)) -\f)dx, 

Jn Jn 
and therefore 

ri-f)= I ($*(-/) + A/)rfx. (34) 
Jn 

The Legendre transform $*(/) of 7r7$(K/(27r)) is given by $*(/) = for 
I/I < 7/2 and $*(/) = 27rfc|/| - 7r7p for k-f - 7/2 < |A;| < k^i + 7/2. 
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Thus flSTj) is equivalent to the problem 



min{^ (|V/|^ + f + 2$*(/) + 2A/) dx; f e H'.iQ)], (35) 

and the limit vorticity is defined in terms of the minimizer / by the formula 

An important role in the analysis of problem f l5^ plays the pointwise 
monotonicity of minimizers in A shown in the following 

Lemma 5. Let fx be the minimizer of l[35\} . then (i) fx < in Q for every 
A > 0, and fa^f/3'i'nflifa> (3. 

Proof. Approximate the function $*(/) by smooth convex functions $J(/), 
where 5 > is a small parameter. Set 



f-s 



Clearly 4>^ G C^(]R) and is convex thanks to the convexity of $*(/)• Let 
the fx be the minimizer of the functional fl35l) with $| in place of $*. It is 
obvious that this minimizer is continuous function. Assume that the function 
f = fi3 — fa has a negative minimum. Subtracting Euler-Lagrange equation 
for from that for we obtain 

-A/ + / + (/3 - a) + (l>|)'(/^) - i^infa) = 

At the minimum point of / we have that —A/ < 0, f < 0, (3 — a < 
and ($|)'(//3) — ($|)'(/a) < 0. Thus we have a contradiction and therefore 
fli^fa- In particular, setting /3 = we get /a < for A > 0. 

The result follows by passing to the limit 5 — )■ 0. □ 



4.1 Weak magnetic fields: zero vorticity 

Let us consider weak magnetic fields h'l^^ = A/e^, such that A > is small. 
It is natural to expect that for such magnetic fields the minimizer fx of the 
problem ( 135|) satisfies —7/2 < /a < therefore $*(/a + v) = in fl for every 
sufficiently small smooth test function, and this implies that fx must solve 
the problem 

'Af = f + Xmn , , 

36 

/ = on 

More precisely, the case of zero vorticity is described by 



13 



Proposition 1. Let fi be the solution of ([23) for A = 1 and 7 is defined in 

m- If 

A < A,,i := ^ (37) 
2 max I ji I 

i/ien the minimizer fx of [3D\] is given by fx = Xfi, and, according to ([23), 
D{x) = 0. 

Proof. Since 

I I (|V/p + /2 + 2$*(/) + 2A/)rfx>i I {\Vf\' + f + 2\f)dx (38) 

for every /, and (13 8 p becomes equality if / is the minimizer of the right hand 
side, / = A/i, the result follows. □ 

4.2 Moderate magnetic fields: simple vortices 

When the parameter A increases A/i is no longer the minimizer of (135|) . the 
latter is rather described by the variational problem 

min m (|V/|2 + f + 47r(|/| - 7/2)+ + 2A/) dx; f G H',{n)j . (39) 



Proposition 2. Let gx be the minimizer of ( T^Pj) and let 

Acr2 := max{A > 0; max \gx\ < 37/2}, (40) 
then the minimizer fx of l[35\} coincides with that of ^3^. 



Proof. The proof is identical to that of Proposition [H Note that we make use 
here of Lemma [5] to get that the minimizer gx of fl39l) satisfies the inequality 
gx > —37/2 in Q when A < Acr2- D 

In order to describe the limiting vorticity function D{x) we have to study 
in more details problem ( 139|1 . 

Proposition 3. If Xai < A < 27r + 7/2 then ([5^ is reduced to the obstacle 
problem 



mm 



^ (I V/r + f + 2A/) dx; f e H'.in), I/I < 7/2} . (41) 



The minimizer fx takes value —7/2 on a set with nonzero measure (2D 
Lebesgue). Moreover, the vorticity D{x) is zero in the domain where f{x) > 
—7/2 and D{x) = (A — 7/2)/(27r) otherwise. 

// max{Acri, 27r + 7/2} < A < Acr2 then D{x) = 1 when f\{x) < —7/2 
and D{x) = if fx{x) > —7/2, where fx is the minimizer of [3^1 . 
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Remark 2. It follows from Proposition\^ that the vorticity D{x) is zero in 
a suhdomain of Q where fx > —7/2 and < D{x) < 1 when fx < —7/2. 
Moreover, there are two scenarios, if Xcri < A < 27r + 7/2 then < D{x) < 1 
in the set where fx{x) < —7/2 (holes with degrees one and zero), and if 
max{Acri, 27r + 7/2} < A < Acr2 then D{x) = 1 when fx{x) < —7/2 (all holes 
in this set have degree one). 

Proof Let Acri < A < 27r + 7/2 then + 47r(|/| - 7/2)+ + 2A/ > 7V4 - 
A7 when / < —7/2. It follows that the minimizer fx of (!39|) satisfies the 
pointwise inequality fx > —7/2. Clearly we also have fx < 0. Thus fx 
minimizes fl44p . If we assume that fx satisfies the strict inequality fx{x) > 
—7/2 for a.e. x E fl, then —Afx{x) + /a(x) + A = for a.e. x E Q 
and therefore fx is the solution of problem fl36|) . Since A > Acri we have a 
contradiction with the pointwise bound fx > —7/2. 

In the case when max{Acri, 7r+7/2} < A < Xct2 we clearly have —Afx{x) + 
fx{x) + A = when -7/2 < /a < 0, and -Afx{x) + fx{x) + X = 2% when 
fx < —7/2. Thus we only need to show that the level set fx = —7/2 has zero 
measure. To this end consider the set W = {x E Q; —7/2 > fx > —7/2 — 5}, 
where 6 > 0. For sufficiently small S the boundary of W can be divided into 
two nonempty parts Si = {x E dW; fx = —7/2} and S2 = {x E dW; fx = 
—7/2 — 6}. Both sets Si and S2 have zero measure. Consider the function U 
such that AU = in the interior of W, U = —7/2 on Si, and U = —7/2 — 6 
on 5*2. By the maximum principle U < —7/2 in the interior of W. On the 
other hand fx'^U (otherwise min{/A, f/} is a minimizer). Thus the level set 
fx = —7/2 coincides with Si and has zero measure. □ 

4.3 Stronger magnetic fields: multiple vortices 

For A > Acr2 vortices with multiplicity two appear. Similarly to the case of 
simple vortices there are two scenarios depending on whether Acr2 < 7/2 + 27r 
or Acr2 > 7/2 + 27r. Define 

Acr3 := max{A > 0; max \gx\ < 57/2}, (42) 

where gx is the minimizer of the problem 

min m (|V/r + f + 47r((|/| - 7/2)+ + (|/| - 87/2)+) + 2A/) dx; 

fEH',{Q)}. (43) 
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Proposition 4. If Xct2 < A < 47r + 87/2 then ( f^^y zs reduced to the obstacle 
problem 

mill |i ^ (I V/p + f + M\f\ - 7/2)+ + 2A/) dx; f G iJo'(^)> I/I < 37/2| • 

(44) 

The minimizer fx takes value —87/2 on a set with nonzero measure (2D 
Lebesgue). Moreover, the vorticity D{x) is zero in the domain where f\{x) > 
-7/2, D(x) = 1 when -87/2 < fx{x) < -7/2 and D{x) = (A - 87/2)/(27r) 
when fx{x) = —87/2. 

// max{Acr2) 47r + 87/2} < A < Acrs then D{x) = when f\{x) > —7/2, 
D{x) = 1 if -87/2 < fx{x) < -7/2 and D{x) = 2 when fx{x) < -87/2, 
where fx is the minimizer of ( [73] ]. 

Remark 3. In the case when Acr2 < A < 47r + 87/2 we see that vortices 
with multiplicities one and two coexist in the set where fx{x) = —87/2, while 
all holes have degrees one in the domain where —87/2 < fx{x) < —7/2 and 
zero degrees in the domain where f\{x) > —7/2. // max{Acr2, 47r + 87/2} < 
A < Acr3 there are three subdomains, where f\{x) > —7/2, —87/2 < fx{x) < 
—7/2 and fxix) < —7/2. All holes in these domains have degrees zero, one 
and two, correspondingly. 

The proof of this result is similar to the previous ones. Further increase 
of the magnetic field leads to vortices with higher multiplicities in nested 
subdomains. 

Proposition 5. There exists a strictly increasing sequence of critical values 
\rj, 3 = 1,2,... such that if X^rj < A < Acr^+i) the limiting vorticity takes 
constant values in subsets Qk \ ^k+i, where Qk = ^k{X), = 0, 1, . . . , j are 
strictly nested sets and Qq = Q. Namely, the vorticity D{x) = in Qq\Qi 
and D{x) = k in flk\ ^^fc+i , k < j — 1. Finally, when x G Qj then there are 
two scenarios: (i) if X < 27rj + (j — 1/2)7 then (j — 1) < D{x) < j otherwise 

(a) D{x) = j. 
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5 Appendix 

The derivation of the lower bound in Section [3] makes use of the concept of 
Young measures. Recall that a Young measure is a parametrized family of 
probability measures nix associated associated to a family of functions 0e(x) 
{(p.-.n^ M) such that 



for every bounded continuous function F{\). It is known (see, e. g., [1], 
[11], [13]) that under some a priory bounds on the sequence of function (/)e(x) 
there exists a family rrix such that fH5l) . 

In this work we consider a sequence of integer valued functions D'^{x). We 
construct a partition of unity /ifc(x) associated this sequence via fl2^ - fl251) . so 
that the corresponding Young measure on 9} is the sum of (5-functions 
centered at integer points, mx{X) = J2 f^k{x)6k{X)- For fixed k the value 
fj,k{x) represents the probability to find a hole Uj with degree dj = kma. 
small vicinity of the point x (i.e. ^k{x) represents the ratio of holes with 
degree in a small neighborhood of x to the total number of holes in this 
neighborhood). 




(45) 
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